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We study the predictive power of Skyrme forces with respect to low lying quadrupole spectra 
along the chains of Sn, Cd, and Te isotopes. Excitation energies and B{E2) values for the lowest 
quadrupole states are computed from a collective Schrodinger equation which as deduced through 
collective path generated by constraint Skyrme-Hartree-Fock (SHF) plus self-consistent cranking 
for the dynamical response. We compare the results from four different Skyrme forces, all treated 
with two different pairing forces (volume versus density-dependent pairing) . The region around the 
neutron shell closure N = 82 is very sensitive to changes in the Skyrme while the mid-shell isotopes 
in the region N < 82 depend mainly on the adjustment of pairing. The neutron rich isotopes are 
most sensitive and depend on both aspects. 

PACS numbers: 21.10Dr; 21.10.Pc; 21.60.Jz 
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I. INTRODUCTION 

A key feature of nuclear excitations are the low-lying 
2 + states. Their properties delivered crucial input for 
developing an understanding of nuclear structure Q, • 
At first glance, they suggest the collective picture of 
the nucleus as a liquid drop which can undergo global 
quadrupole oscillations and which freeze under certain 
conditions to a stable rotator. This view has been for- 
mulated in terms of the Bohr-Hamiltonian which estab- 
lishes a collective dynamics in the five quadrupole degrees 
of freedom |3( . The parameters of the collective Hamilto- 
nian have to be adjusted phenomenologically, see e.g. the 
applications in 0. The collective approach has been re- 
vived with the interacting boson model (IBM) which has 
found widespread application and proven to be extremely 
useful in sorting nuclear low-energy spectra [f| . 

The collective picture is seemingly in contrast to the 
microscopic view which sees the nucleus as consisting out 
of shells of single nucleons arranging themselves in a com- 
mon mean field 0, Q • The views can be unified by the 
concept of a deformed mean field which establishes a re- 
lation between single-particle shell structure and global 
deformations |g, |9( ■ The collective motion is then under- 
stood as vibration (or rotation) of the mean field similar 
to the Born-Oppenheimer method for describing molec- 
ular vibrations. The connection is established on a for- 
mally sound level by the generator coordinate method 
(GCM) 0, ^3 which describes collective dynamics as 
coherent superposition of a continuous set of deformed 
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mean-field states, called the collective path. The GCM 
within the Gaussian-Overlap-Approximation (GOA) al- 
lows to establish contact between the microscopic foun- 
dation and a collective Bohr-Hamiltonian (12L Il3| . Start- 
ing from the GCM, the lines of applications spread enor- 
mously. There are, on the one hand, fully fledged GCM 
calculations which skip the collective Hamiltonian as in- 
termediate level and compute low-energy spectra directly 
from the coherent superposition of the collective path; 
these sophisticated calculations imply exact projection 
for the conserved quantities, as particle number, angular 
momentum, and center of mass; there are many pub- 
lished results around, we mention here 0, as two 
recent examples. On the other hand, one finds several 
approximations to the microscopic computation of the 
Bohr-Hamiltonian; most applications hitherto employ a 
phenomenological shell model to describe the deformed 
mean field, see e.g. ^(J. There are also several self- 
consistent calculations along that line, for an early ex- 
ample see |T3 and for more recent achievements [lSI Il9j . 

An alternative direction of development remains at the 
microscopic mean-field description and makes it man- 
ageable by restricting considerations to small amplitude 
motion. This yields the much celebrated random-phase 
approximation (RPA) which has its stronghold in the 
description of giant resonances, see e.g. The ap- 

propriate extension to non-closed shells with pairing is 
the quasi-particle RPA (QRPA) which has only recently 
been developed up to a rigorously self-consistent level 
[2l|. The QRPA describes formally the whole excitation 
spectrum, including the low lying 2 + states, and it opti- 
mizes all states automatically. It assumes, however, small 
amplitudes, i.e. harmonic motion. This is perhaps legit- 
imate close to magic nuclei but somewhat dubious else- 
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where. The above mentioned theories for large amplitude 
collective motion concentrate on the lowest state only 
but try to take into account all effects of anharmonic- 
ity due to soft potential energy landscapes and shape 
isomerism. The fully fledged adiabatic time-dependent 
Hartree-Fock method (ATDHF) provides an unambigu- 
ous optimization scheme for the large-amplitude collec- 
tive path [l2j,|22j. However, that rather involved scheme 
has not yet been used for heavy nuclei as we are going to 
study here. We use presently ATDHF only to compute 
the self-consistent collective mass and employ the more 
intuitive constraint Hartree-Fock to generate the path. 

The connection from a microscopic Hamiltonian to col- 
lective spectra via a large-amplitude collective path is 
well established by virtue of the GCM. An open problem 
is the microscopic input. Self-consistent nuclear mean- 
field models employ effective energy functionals as e.g. 
the Skyrme-Hartree-Fock method, the Gogny force, or 
relativistic mean field, for a recent review see |23j. These 
are empirically adjusted to nuclear ground state prop- 
erties of stable nuclei. There exists a large manyfold 
of equivalent parameterizations which provide compara- 
ble ground state properties but can differ substantially 
in p redictions to exotic nuclei or resonance excitations 
[23L |24| . It is by no means guaranteed that all mean-field 
parameterization produce at once the correct collective 
low-energy vibrations. The contrary is to be expected, 
namely a broad span of predictions amongst which only 
a few parameterizations deliver a satisfying spectrum. 
To phrase that positively: low-energy vibrations provide 
useful information for a better selection of mean-field pa- 
rameterizations. We aim here at a first exploration of 
the connection between mean-field parameterizations and 
emerging low-energy spectra. We do that for the Skyrme- 
Hartree-Fock approximation by comparing the results of 
several different Skyrme forces and pairing recipes. 

It is obvious that such systematic studies need to con- 
fine the subject and the method in order to keep things 
manageable. As test cases, we consider the lowest 2 + 
state in the chain of Sn isotopes and its even neighbors 
Cd and Te. These share basically one type of collec- 
tive motion being predominantly soft vibrators. For the 
practical technique, we employ GCM-GOA through a mi- 
croscopically computed Bohr-Hamiltonian. For reasons 
of simplicity, the microscopic information is computed 
along axially symmetric shapes and interpolated into the 
full space of quadrupole degrees of freedom. This ap- 
proximation allows large scale scans and it is acceptable 
for soft vibrators as they are studied here. 



II. FORMAL FRAMEWORK 

A. Underlying microscopic model - input 
parameters 

As starting point, we take a microscopic mean-field 
theory at the level of the Skyrme-Hartree-Fock model 



augmented by pairing in BCS approximation plus Lipkin- 
Nogami correction for approximate particle number pro- 
jection. This is a standard approach in nuclear structure 
physics. We refer the reader to |2^| for a detailed descrip- 
tion of the energy functional and subsequent mean-field 
equations. We recapitulate here only briefly the spec- 
trum of variants of that model which will play a major 
role in the following discussions. 

The mean-field part is determined by the Skyrme en- 
ergy functional E^y( y p v ^T v ^J v ,j l j,a v ) which depends on 
the local density p v , kinetic-energy density t„ , spin-orbit 
density J v , current j„ and spin-density a v and where q 
means protons or neutrons. The functional form is basi- 
cally settled since two decades [2j| with minor extensions 
in later stages (e.g. ptl 127^). However, there exists a 
great variety of actual parameterizations for the Skyrme 
energy functional. Most of them provide a high-quality 
description of nuclear bulk properties as binding ener- 
gies and radii. They differ in details as, e.g., isovector 
forces or surface properties. We are going here to ap- 
ply the Skyrme functionals to a regime far from what 
had been considered in the fits. It is thus important to 
explore a minimal variation of parameterizations within 
the Skyrme framework. We will consider here: SkM* as a 
widely used traditional standard [25| , Sly6 as a recent fit 
which includes information on isotopic trends and neu- 
tron matter |28( , SkI3 as a fit which maps the relativistic 
isovector structure of the spin-orbit force and takes care 
of the surface thickness [2(|, and SkO [2^| as a recent 
fit relying an the same fit data as SkI3 but with addi- 
tional constraint on the two-nucleon separation energies 
around 208 Pb and with a better adjusted asymmetry en- 
ergy. That selection contains a large span of effective 
masses: SkI3 <-» m*/m — 0.6, SLy6 «-* m*/m = 0.7, 
SkM* <-> m*/m = 0.8, and SkO <-► m*/m = 0.9. The 
effective mass has an influence on the level density near 
the Fermi surface which, in turn, may have an effect on 
the low-energy collective states. There is also a differ- 
ence in the isovector and spin-orbit properties. Besides 
the effective mass and asymmetry, the bulk parameters 
(equilibrium energy and density, as well as incompress- 
ibility) are comparable. 

The second key ingredient is pairing. A present-day 
standard is to use a zero range pairing force often called 
volume pairing. We will use the notion (S-interaction (DI) 
pairing. A widely used variant for the pairin g fo rce is a 
density-dependent delta-interaction (DDDI) 30] . Both 
recipes are summarized as 



vS m) 8(n-r 2 ) 



V,) 



(DDDI) 



S(n - r 2 ) [1 - p(r)/po] 



(1) 



The pairing strengths or yJ DDDI ^ are adjusted to odd- 
even staggering of binding energies in a few representa- 
tive semi-magic nuclei (Sn and Pb isotopes, N=82 iso- 
tones). The adjustment is done for each force separately 
because the much different effective masses call for dif- 
ferent pairing strengths. The actual values used here are 
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SkM* 
SLy6 
SkI3 
SkO 


279.1 259.0 
298.8 288.5 
335.4 331.6 
253.0 269.0 


990.0 802.0 
1053.1 864.2 
1233.0 996.0 
1007.4 893.7 



TABLE I: Pairing strengths for the two pairing recipes and 
for the Skyrme forces used in this paper. The strengths are 
given in units of fm— 3. 

given in table D The pairing recipe is to be augmented 
by a cutoff in single particle space. We use a smooth 
cutoff with a Woods-Saxon profile in the single particle 
energies. The switching energy is chosen such that the 
pairing space covers 1.6iV 2 / 3 particles above the Fermi 
energy, for details see In order to explore the influ- 
ence of the pairing recipe, we will also discuss deliberate 
rescaling of the pairing strengths. 

B. Deduced collective dynamics 

The mapping from the microscopic to a collective 
description is performed with the generator-coordinate 
method (GCM). This is a much celebrated method in nu- 
clear structure phy sics, for a review see e.g. an d for a 
brief summary |23]. We outline here the basic steps and 
provide a more detailed compact account in appendix lAl 

The stationary mean field equations as such provide 
only a few well isolated states, preferably the ground 
state and perhaps some isomers. Each state is char- 
acterized by one BCS wavefunction |$) which is com- 
posed of a set of single-nucleon wavefunctions together 
with their occupation amplitudes. In order to describe 
motion, one needs to consider a time-dependent mean- 
field theory, in the nuclear community often called time- 
dependent Hartree-Fock (TDHF). Large amplitude col- 
lective motion is related to low energy excitations, thus 
slow motion. This justifies the adiabatic limit known as It 
yields at the end, a collective path where q stands 

for continuous series of deformations, predominantly of 
quadrupole type because nuclei are softest in that degree 
of freedom. The dynamical aspect is added in first order 
of collective velocity, i.e. in terms of linear response to 
a collective displacement. It is a widely used approxi- 
mation to determine the collective path from quadrupole 
constrained Hartree-Fock-BCS (CHF). 

The link between microscopic description and collec- 
tive dynamics is established by the collective path 
where q stands for a continuous series of deformations, 
predominantly of quadrupole type because nuclei are 
softest in that degree of freedom. A systematic theory 
for an optimized collective path is provided by adiabatic 
TDHF (ATDHF) H3, H3, |3J| The overwhelming major- 
ity of practical applications simplifies the construction 
by using a simple quadrupole constraint Hartree-Fock to 
produce the The path, once established, serves as a 
basis along which the collective motion expands. The cor- 



responding microscopic state is described as a collective 
superposition |^) = J dq\$ q )f(q). The state is opti- 
mized by a variational principle. This is the fully fledged 
GCM which can be those days attacked in a straightfor- 
ward numerical manner, see e.g. |l4l Il5l l35| . However, 
that is still a very demanding task and not so well suited 
for broad surveys as we intend it here. We use here as a 
simple, efficient, and reliable shortcut the Gaussian over- 
lap approximation (GOA) which parameterizes the norm 
and Hamiltonian overlap in terms of Gaussians, e.g. for 
the norm overlap as = exp (-X(q - q') 2 /4). It 

provides an acceptable approximation, particularly for 
medium and heavy nuclei The GCM-GOA yields 

at the end a fairly simple collective Hamiltonian where 
the collective potentials and masses are unambiguously 
computed from the microscopic energy functional and 
the collective path. Quadrupole motion has five degrees- 
of- freedom 0. [!■ 0]. The emerging collective Hamilto- 
nian thus has the form of a generalized Bohr-Hamiltonian 
while its potentials and masses are computed from mi- 
croscopic input 0, 0, . 

The practice of GCM-GOA is a bit involved, see the 
appendix for a few more details. We summarize here 
the steps: The energy expectation value along the path 
yields a raw collective potential V(q). The collective 
mass and moments of inertia are obtained by dynami- 
cal linear response about a given point at the path often 
called self-consistent (or ATDHF) cranking [3^| ; as an ap- 
proximation, Inglis cranking is used whenever justified. 
Zero-point energy corrections to the potential are com- 
puted from these masses and the collective fluctuations 
(quadrupole, angular momentum) of the states In 
fact, a topologically corrected GOA is used to allow a 
numerical robust computation of potentials and masses 
in the intrinsic frame (defined by a diagonal inertia ten- 
sor) |36t l37| . This provides an interpolation scheme to 
connect safely the near spherical shapes with larger de- 
formations. 

After all, we restrict the microscopic calculations to 
axial symmetry. The fully five-dimensional quadrupole 
dynamics is recovered by interpolation of the collective 
potential and masses between prolate and oblate shapes 
into the triaxial plane. This approximation saves two 
order of magnitude computation time and thus allows 
the large scale systematics as we intend it here. On the 
other hand, it is well justified at and in the vicinity of 
the spherical shape. The test cases for the present study 
are Sn isotopes and its even-even neighbors Cd and Te 
which are predominantly soft vibrators around spherical 
mean. Moreover, we confine the study to the first excited 
2+ state (and occasionally to the + ground state) which 
both are not very sensitive to details in the triaxial plane. 
All that considered, the triaxial interpolation is a useful 
and legitimate approximation for the intended systematic 
explorations. 
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FIG. 1: The raw collective potential V and the effective po- 
tential V including the zero-point energy corrections \ A6I . 
both drawn as function of the intrinsic axial quadrupole mo- 
mentum 020- Test case is 132 Sn computed with DI pairing 
and SkI3. The position of the ground state and the first 
2 + state are indicated by horizontal bars. The difference be- 
tween the minimum of V and the 0+ energy is the correlation 
energy AE COII . 



C. An example for potentials and masses 

In order to exemplify details of the calculations, figure 
^ shows for the case of 132 Sn the collective potential be- 
fore and after zero-point energy correction (ZPE) 1AGII . 
The ZPE induce obviously a strong global down-shift in 
energy because the spurious energy content from collec- 
tive fluctuations in the is subtracted. Moreover, they 
may change the shape of the potential. The corrected po- 
tential has its minimum at a slightly deformed position 
although the doubly magic 132 Sn is a perfectly spherical 
nucleus in a pure mean-field description (see the well de- 
fined spherical minimum in V). This is the same effect 
as happens in variation after rotational projection (for a 
model discussion see |36(): knowing that the projection 
restores spherical shape anyway, the system takes advan- 
tage of a small deformation to acquire correlation energy. 
It is comforting that we see the same effect here because 
our treatment of quantum correction should include a 
good approximation to rotational projection. 

Figure ^ also indicates the position of the 0^ ground 
state and the first excited 2 + state. The lies above 
the bottom of the intrinsic potential V as it should be, 
to account for the correct physical zero-point energies, 
but it stays below the minimum of the raw potential V 
because the larger spurious zero-point energy had been 
subtracted before adding the physical one. The net ef- 
fect is a correlation energy AE COTI which expresses that 
the collectively correlated ground state is better bound 
than the mean-field ground state. The 2 + state lies, of 
cause, above the Oq state. The quantity of interest here 
is the excitation energy E(2 + ) which is computed as the 
difference between the total 2 + energy and the Oq energy. 
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FIG. 2: Energies E(2 + ) and B(E2) | values (= 
|(0 + |Q2m|2 + )| 2 ) along the chain of Sn isotopes calculated us- 
ing the four different Skyrme interactions as indicated. The 
experimental results are taken from [3S|. 



force 


SkI3 SLy6 SkM* SkO 


E{2 + ) [MeV] 

spectral gap protons [MeV] 
spectral gap neutrons [MeV] 
m* jm 


4.36 3.76 3.94 2.41 
6.6 6.2 6.4 6.0 
6.6 6.0 5.4 4.0 
0.6 0.7 0.8 0.9 



TABLE II: Comparison of the E{2 + ) energies in 132 Sn for 
the various Skyrme forces with the spectral gaps in 132 Sn and 
the effective mass m*/m associated with the forces. 



III. RESULTS 
A. Results for the chain of Sn isotopes 

1. Variation of forces 

Figure [5] shows the 2 + excitation energies and tran- 
sition strengths along the chain of Sn isotopes for the 
four chosen SHF parametcrizations and in comparison to 
experimental data. At first glance, we see that all cal- 
culations hit the right order of magnitude. They also 
reproduce the increase of the E 2 + at the shell closure 
N = 82. At closer inspection, however, we see interest- 
ing differences and mismatches in detail. 

Let us first concentrate on the doubly magic case of 
N = 82. Shell effects directly related to the SHF forces 
should dominate here. And indeed, we see that the E(2 + ) 
are closely related to spectral properties. Table ITT1 shows 
the E(2 + ) energies in comparison to the spectral gaps 
of protons and neutrons in 132 Sn. The spectral gaps are 
the energy difference between the highest occupied single 
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particle orbital and the lowest unoccupied orbital (known 
as HOMO-LUMO gap in molecular physics). The neu- 
trons show always the smaller gaps and these lowest 1- 
particle-l-hole (lph) transitions take the lead in the com- 
position of the lowest 2+ state. Correspondingly, both 
quantities share the same trends. The spectral gap, in 
turn, is related to the effective mass m* /m of the forces. 
We see that also in table |H] where low m*/m correlate 
to large gaps and vice versa. But the step down to the 
rather low spectral gap for SkO is much larger than the 
step up in effective mass. Here we see also an interfer- 
ence from the very strong isovector spin-orbit force of 
SkO, another important contributor to shell effects. The 
tablelTll furthermore, demonstrate the effect of the resid- 
ual interaction in that the E(2 + ) are generally 1.5 MeV 
below the lowest lph energy. This allows to postulate 
a simple criterion for the selection of forces: the lowest 
spectral gap in 132 Sn (and other doubly magic nuclei) 
should stay safely above the experimental E(2 + ), which 
is 4.04 MeV for 132 Sn. The force SkO clearly fails in that 
respect. The reason is that SkO was fitted to match the 
two-nucleon shell gaps at doubly magic 208 Pb already at 
the level of pure mean- field calculations |2^. Meanwhile, 
it has been shown that collective correlations reduce the 
two-nucleon shell gaps by 1-2 MeV [3?! . The fitting strat- 
egy of SkO thus squeezes the spectral gap too much with 
the obvious consequence that the collective spectra are 
spoiled throughout. This mismatch is thus a strong hint 
on the inner coherence of Skyrme forces connecting the 
various observables. 

Far away from the magic N = 82, one expects that the 
pairing gap dominates the E(2 + ) energy. The pairing 
force was tuned in the same way to the odd-even stag- 
gering in Sn isotopes and iV = 82 isotones. Thus the 
pairing gap is about the same for all four Skyrme forces 
in the well pairing region [N < 80) . We see indeed com- 
parable energies for the three forces, SkM*, SLy6, and 
SkI3, which also hit very nicely the experimental values. 
The force SkO, however, produces systematically larger 
energies out there. This shows that shell effects (here 
probably from the spin-orbit force) have also some in- 
fluence. Different relations are seen in the other pairing 
regime for the neutron rich, exotic nuclei above N = 82 
where SkM* shows always the largest energies. This is 
at the same time a region of weak binding. This causes 
a strong interplay of shell effects and pairing which are 
not easily disentangled. 

Large differences are seen in the immediate vicinity 
of N = 82. All forces, expect SkO, reproduce nicely 
the sudden step down from N = 82 and the asymmetry 
around N = 82, namely the fact that N = 84 has lower 
E{2+) than N = 80. But the results differ in the trends 
for N — 78. The case SkI3 follows nicely the smooth 
experimental trend while SkM* and SLy6 show a spike. 
That is compensated at N = 70 where now SkI3 has a 
spike. In all these cases, we found that a larger E(2 + ) is 
related to a somewhat lower neutron level density at the 
Fermi surface. 



A much more critical test than excitation energies are 
the associated transition probabilities, the B(E2) values. 
One is usually happy to describe them within a factor 
of 2 or so, and often the concept of effective charges is 
introduced to achieve a fine tuning |40j| . The lower panel 
of figure |21 shows the B(E2) J \ values for the transitions. 
For N < 82, they are similar for all four forces in spite 
of the sometimes very different energies. But they all 
differ from the experimental data by about a factor of 
two. The positive aspect is that the theoretical results 
come so close at all in view of the fact that the transi- 
tion strengths are always much more demanding. The 
remaining mismatch can have various origins: 1) We use 
simply a quadrupole constraint to generate the collec- 
tive path instead of the variationally optimized ATDHF 
prescription [3j, |33, |34( . 2) We use the raw quadrupole 
expectation value rather than the fully mapped collective 
image (see appendix IA 6(1 . 3) The effective energy func- 
tional is not fully suited to compute transition moments 
and effective charges had to be added for a correct de- 
scription jijj • Which one of these approximations is most 
responsible, has yet to be explored. Anyway, the results 
are not untypically bad because almost all microscopic 
approaches have a hard time with an exact reproduction 
of transition moments. 

The B(E2)1 in the region around the doubly magic 
132 Sn surprisingly shows basically no differences between 
the forces, just in a region where the energies differ most. 
A very interesting point is the magic 132 Sn. Naive mod- 
els predict a dramatic drop in the B(E2) at the magic 
point. A model study taking care of the residual inter- 
action and cross-talk between the neutron and proton 
quadrupole vibrations predicts that the B(E2) should, 
quite oppositely, have a peak at 132 Sn 41]. Our cal- 
culations confirm these estimates at a qualitative level, 
namely to the extent that we also do not find any deep 
dip in the B(E2). In our cases, the residual interaction 
was obviously not large enough to turn that into a peak. 
But these are quantitative details of the employed forces. 

The largest differences between forces for the B(E2)^ 
values are seen in the deep exotic regime A > 132. We 
are sure that information about B(E2) in that region 
would be valuable. But before one can exploit that, one 
has to understand (and possibly remove) the systematic 
ovcrcstimation still seen on the low- A side. 



2. Variation of pairing recipes 

Figure[3]shows the collective spectra along the Sn chain 
for SkI3 computed with different pairing prescriptions. 
We have added in the lowest panel some information 
about the internal pairing structure, namely the aver- 
age neutron-pairing gaps A = J2 a u a v a A aa / J2 a u a v a 
which are deduced from spectral properties of the given 
nucleus at th spherical shape (usually the minimum n the 
PES) and which, nonetheless, provide a simple measure 
for the pairing gap deduced from odd-even staggering 
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FIG. 3: Energies E(2+), B{E2)] values (= |{0 + |Q 2 m|2+)| 2 ), 
and average neutron-pairing gaps at spherical shape along the 
chain of Sn isotopes calculated with SkI3 and three different 
pairing recipes: standard DI (circles), DDDI (boxes), and DI 
with 25% enhanced strength. The experimental results are 
taken from |38|. 



|3lj . The deliberately changed pairing strengths (boxes 
versus circles for enhanced pairing and rhombus versus 
circle for reduced pairing) has an obvious effect. The 
pairing gap is increased or reduced and subsequently, the 
E(2 + ) energies change in the same direction. The effect is 
most pronounced in the regions sufficiently far off N = 82 
where we expect a dominance of pairing in the collective 
spectra. The relative changes in excitation energies and 
pairing gaps are much larger than the change in pairing 
strength. Moreover, the excitation energies behave gen- 
erally similar to the pairing gaps. This demonstrates that 
the low-energy spectra in soft vibrators provide valuable 
information about the pairing strengths. On is tempted 
to use that for an immediate tuning of the strengths. We 
run however, in some conflict, because E(2 + ) energies 
and B(E2)1 do no coincide at the same strength. More- 
over, one has to keep in mind that the information from 
lying states is still mixed with effects of the mean field. 
This is a general feature of nuclear structure, even for 
the odd-even staggering which is usually taken as bench- 
mark for pairing properties [42| . The interference of shell 
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FIG. 4: As figure |21 but now with DDDI pairing. 



effects can be seen here particularly well from some ir- 
regularity at N = 68. The reduced pairing produces for 
N = 68 a strongly deformed ground state which results 
in a sudden drop of the E(2 + ) energy accompanied by 
a strong peak in B(E2) J \. The average gap shows no 
dramatic reaction because it remains related to the now 
irrelevant spherical shape. The DDDI pairing stays in 
most cases more or less close to the results of DI pairing 
which is somewhat expected because it is tuned to the 
same average pairing gap. However, DDDI pairing reacts 
differently to shell effects as can be seen in the vicinity 
of the shell closure N = 82 and for the weak shell closure 
which SkI3 produces at N — 72. The DI pairing seems to 
comply better with data. But that holds in connection 
with the particular shell structure of SkI3. Much more 
systematic investigations with varied forces and in other 
region of the nuclear chart are necessary before drawing 
any conclusion like that. The B(E2)1 values shown in the 
lower panel of figure [3] show generally smooth trends, ex- 
cept for N > 86 where an increasing trend sets on which 
is related to the increasing softness of these neutron rich 
isotopes. The sensitivity to varying pairing recipes is sim- 
ilar to what we have seen when varying the forces: They 
vary little for N < 82 and more significant differences 
appear in the far exotic regime N > 82. 

For completeness, it is worth to look at the perfor- 
mance of DDDI pairing also for the other Skyrme forces 
in the survey. That is shown in figure 0] It has to be 
compared with figure [21 Similarities and differences are 
about the same for all shown forces. The average exci- 
tations in the well pairing regime are comparable. The 
small fluctuations about the average trends appear for 
DDDI at different places than for DI pairing. The most 
pronounced difference to DI is seen for the E{2 + ) energy 
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next to the magic neutron number, i.e. for TV = 80 and 
84. DI pairing reproduces the steep experimental drop 
while the DDDI results make all a somewhat less dra- 
matic step. The B(E2)1 values shown in the lower panel 
are very similar to those from DI pairing. They seem to 
be here the more robust signal. We conclude from this 
results that the B{E2) J \ are insensitive to pairing while 
they are the much more sensitive observable in other re- 
spects, e.g. in its dependence on the force. 

A hnal comment on the B(E2) J \ values: Good vibrators 
and well-developed rotators are distinguished by the fact 
that the lowest 2+ state exhausts the quadrupole sum 
rule in collective space. The test for this feature is the 
comparison of the variance (= quadrupole sum rule) with 
the B{E2) value 



<0+|Ql o |0+> 



EK 2 n IQaolO" 1 



(2) 
the 
the 



where (...) means the average in collective space, CP 
ground state, 2+ the spectrum of 2 + states, and 2^ 
lowest 2 + state. We have checked that and found that 
there is generally good exhaustion of the variance by the 
lowest 2 + . The collective potential V(a2o) of all these 
isotopes excludes rotors, so that we can conclude that 
they are good vibrators. 



3. Effect particle-number correction 
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FIG. 5: Energies E(2+) and B(E2) | values (= 
|(0 + |Q2m|2 + }| 2 ) along the chain of Sn isotopes computed 
with SkI3 and DI pairing, once with (solid line, full circles) 
and once without (dashed line, open circle) particle number 
restoration as outlined in section |A 41 The Experimental re- 
sults (dotted lines, full squares) are taken from pjgll . 



The collective Schrddinger equation contains the 
particle-number correction with iVcoii as discussed in sec- 
tion IA 41 It is interesting to check the impact on col- 
lective properties. This is done in hgure There is 
minimal difference for the B{E2) A \ values. The main ef- 
fect is seen for energies in the region of weak pairing, i.e. 
at and around shell closure. In fact, the particle-number 
corrected treatment seems to me a bit more sensitive to 
shell structure as can be seen from the fact that calcu- 
lations without the correction shows generally smoother 
trends. But this statement should be taken with a grain 
of salt. The differences are anyway not very dramatic 
in view of the effects we see when comparing forces and 
pairing recipes. 



Skl3(DI) 
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B. Results for the isotopes of Cd and Te 

The Sn isotopes have a magic proton number Z = 50. It 
is interesting to have a look at its even neighbors, Cd with 
Z = 48 and Te with Z = 52. As a first impression, we show 
in figure[S]a direct comparison for a few selected isotones. 
The effect is obvious. The step from the magic proton 
number to the non-magic ones reduces once more the 
E{2 + ) energies by a substantial factor, fully in agreement 
with the experimental findings. The quadrupole mode in 
Cd and Te is much softer than in Sn where the magic 



FIG. 6: Systematics of the energies of the hrst excited 2 + 
states calculated with the microscopic Bohr-Hamiltonian 1A9II 
using the interaction SkI3 for the nearest even-even neighbors 
of 116 Sn. The Experimental results are taken from |3a| . 



proton number enhances the rigidity of the whole mode 
due to a strong residual proton-neutron force 0] . 

A summary of results for the Cd and Te isotopes using 
the four different Skyrme interactions SkI3, SLy6, SkM*, 
and SkO is displayed in Fig. □ For Te, the E{2+) from 
different forces are very close to each other and to the 
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FIG. 7: Energies E{2 + ) and B(E2) | values (= 
|(0 + |Q2m|2 + )| 2 ) along the chain of Cd isotopes (upper pan- 
els) and Te isotopes (lower panels) calculated with different 
Skyrme forces as indicated. The experimental results are 
taken from |3Sl |. 



experimental data up to N = 76 and show again larger 
differences near the shell closure TV = 82. The case for 
Cd is similar showing, however, an earlier onset of dif- 
ferences. It is noteworthy that the results for SkO reside 
well amongst the other forces, although it behaved dra- 
matically different for the Sn isotopes. This is related 
to the fact that shell effects are somewhat suppressed in 
Cd and Te because these have non-magic proton num- 
ber. Figure [7] shows also as complementing information 
the B(E2) values for the Cd and Te isotopes. The differ- 
ences between the theoretical and the experimental val- 



FIG. 8: As figure Q but now for DDDI pairing. 



ues are for Cd in the same order of magnitude like in 
the case for the Sn chain. They tend to be much less for 
the Te isotopes for N < 78 where deviations stay below 
20%. This looks like a remarkable agreement. But both, 
agreement for Te and disagreement for Cd and Sn have 
yet to be understood in detail. 

We have seen in the Sn chain that the step from DI to 
DDDI pairing makes most differences next to the neutron 
shell closure. One has to suspect that a similar feature 
appears next to the proton shell closure. The Sn chain 
reside at N — 50 which is a closed proton shell. Thus the 
neighbor chains for Cd an Tc are in the most sensitive 
regime and we expect visible differences. The results on 
the low lying 2 + states along Cd and Te are shown in 
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figure|Hl We see indeed that the E(2 + ) energies are larger 
with DDDI, particularly near the neutron shell closure at 
TV = 82. The differences become again negligible far out 
in the well pairing regime. And, as for Sn, the B(E2) 
values are totally insensitive. It is also clear that the 
DI results, here and in the Sn chain, are closer to the 
data. The same effect was already seen for the neutron 
channel in figure [21 The step of the E(2 + ) when moving 
away from a magic number is softer for DDDI than for 
DI. This produces somewhat to high E{2 + ) for N = 80 
and 84 in figure [3] and here for Z — 48 and 52 in figure 
|7| vs. figure |H1 We just remark this observation. It is to 
early to draw far reaching conclusions on the validity of 
DI versus DDDI. The difference is seen in the worst case, 
namely a nucleus in the weak pairing regime where we 
are not yet sure that the present pairing treatment (BCS 
plus Lipkin-Nogami) is fully appropriate. 

C. The isotope shifts for the Sn isotopes 

The ground state solution of the collective Schrodinger 
equation provides the collective ground state correla- 
tions. One can compute the correlation effect on any 
one-body observable with the help of the collective map 
as outlined in section IA 61 Figure shows results for 
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FIG. 9: Isotope shifts of the charge r.m.s radii, S(r 2 m!i ) 66 ' N , 
relative to lle Sn. Upper panel: comparison of results from 
the four forces SkI3, SLy6, SkM*, and SkO including collec- 
tive ground state correlations. The expectation value (r 2 ) is 
calculated here according to lA 61 Lower panel: Comparison 
of pure mean field result with those including correlations for 
the force SkI3. The experimental data is taken from 

the systematics of charge r.m.s. radii drawn in terms of 



isotope shifts relative to 116 Sn. The r.m.s. radii shown 
in the upper panel here are taken from the correlated 
ground state which includes the collective shape fluctu- 
ations. At first glance, all forces reproduce the experi- 
mental trend very well where data are available. There 
remain small but significant differences between the four 
forces in that range amongst which SkI3 comes generally 
closest to the data. The similarity of the trends persists 
to isotopes with larger neutron numbers. Substantial dif- 
ferences develop at the upper end for the shown chain, 
i.e. for N > 86. Not surprisingly, this is the regime of 
exotic nuclei because the generally soft binding amplifies 
small differences in shell structure. It is surprising, how- 
ever, that these differences develop so late. The regime 
of similarity reaches well beyond the N = 82 shell clo- 
sure. This is due to the smoothing features of the shape 
fluctuations. The lower panel demonstrates the effect of 
ground state correlations for that observable. There are 
practically no visible effects as compared to pure mean 
field calculations, as one could have expected for such a 
chain of semi-magic nuclei [2(| . Correlation effects be- 
come visible again at the upper edge of the chain where 
the general softness of the deeply exotic nuclei does also 
allow for larger shape fluctuations. 



IV. CONCLUSION 

We have investigated the predictive power of nuclear 
effective forces for describing low-lying collective states 
considering as test cases the chain of Sn isotopes as well 
as its even neighbors Cd and Te. As a particular example 
of such an effective force we used the Skyrme-Hartree- 
Fock scheme augmented by a short-range pairing force. 
To that end, we used a representative sample of differ- 
ent Skyrme forces as well as two different pairing models 
(volume pairing versus density-dependent pairing). The 
spectra of collective quadrupole vibrations were com- 
puted in a two step procedure: First, mean-field cal- 
culations with quadrupole constraint and self-consistent 
cranking where performed which provide the microscopic 
input for a collective Hamiltonian in terms of potentials, 
masses and moments of inertia. Second, the collective 
Schrodinger equation thus obtained is solved in the space 
of the five quadrupole coordinates. Care has been taken 
to subtract correctly the zero-point energies from spuri- 
ous collective fluctuations in the mean-field states and to 
respect the topology of the quadrupole space. As a sim- 
plification, we use axially symmetric mean-field calcula- 
tions and interpolate triaxial properties between prolate 
and oblate shapes. This is an acceptable approximation 
for the nearly spherical soft vibrators considered in the 
present survey. Pairing is treated at the level of BCS with 
particle number correction in terms of the Lipkin-Nogami 
approximation. A final fine-tuning of the average particle 
number is performed also for the collective states. 

We find three regimes: collective properties are domi- 
nated by the pairing gap for N<76, they are dominated 
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by the spectral gap of the neutron level for 76 < N<82, 
and a sudden transition to prolate ground state defor- 
mations emerges for N > 82. In the pairing dominated 
regime, the results for the 2 + excitation energies depend 
mostly on the pairing strength and only weakly on the 
Skyrme forces (with exception of the force SkO which 
behaves a bit strange in consequence of the constraint on 
two-nucleon shell gas in the fit). As all pairing models 
used here were fitted to the odd-even staggering in the 
Sn region, we find generally nice agreement with exper- 
imental data in that mid shell region. In the shell-gap 
dominated regime, on the other hand, a strong sensitiv- 
ity to the Skyrme force develops due to a strong relation 
to the spectral gap which, in turn, depends sensitively on 
the effective mass. These features persist into the regime 
above N = 82. The result for the transition probabili- 
ties (the B(E2) values) show larger deviations from the 
data (up to a factor of two). This is not surprising be- 
cause B(E2) values are generally more demanding to any 
model. 

To summarize, we have shown that Skyrme forces have, 
in principle, the capability to describe low lying collective 
spectra. In practice, the success depends on the actual 
parameterization used. Turning the argument around, 
we find that a systematic investigation of collective spec- 
tra delivers extremely useful information for the selection 
of parameterizations and the development of improved 
effective forces. This calls for more systematic investiga- 
tions. 
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APPENDIX A: MICROSCOPIC COMPUTATION 
OF COLLECTIVE OPERATORS 

1. The deformed mean field 

Microscopic basis is a self-consistent mean-field accord- 
ing to SHF with DI or DDDI pairing (for details see p^l. 
The SHF-BCS equations describe the nuclear state in 
terms of a set of single particle states ip n with associated 
BCS occupation amplitudes v n . These together compose 
the BCS state |i>) = ]J n (u n + v n a+ a~^ n ) |vac) where 

u n = v% ■ As a synonym for its content we denote 

it by |$) = {(p n ,v n }. In practice, we go somewhat be- 
yond the BCS scheme by using the Lipkin-Nogami (LN) 
approximation for particle number projection ^5[. The 
mean-field equations can be summarized as 



(Ala) 



Ct2r 



?20 = r 2 y 2 o/cut(r) , 

4tt ($ Q20 |r 2 F 2 m|$o 



At 



(Alb) 
(Ale) 



The h is a two-quasiparticle operator which itself depends 
on the state on which is acts. The actual form is obtained 
by a functional derivative of the given energy functional. 
The N v is the operator of proton- or neutron- number. 
The Fermi energy e-p^ is chosen to tune the correct par- 
ticle number in the average, i.e. 



£F,, 



($\N U \$)=N V 



(Aid) 



For simplicity we write in the following one particle- 
number term as representative of both. The term oc 
accounts for the approximate particle-number projection 
and its parameter £2,1/ is given according to the LN recipe 
taking properly into account the feedback from the mean 
field to the variances |4^. The LN scheme performs also 
an approximate variation after projection. This yields 
a finite pairing gap under any conditions, even at shell 
closures. And this is the feature what we need to have 
a smooth evolution of the gap along the collective defor- 
mation path. Pure BCS can lead to discontinuities which 
lead to discontinuities in the collective Hamiltonian. 

The stationary mean field equations without constraint 
provide only a few well isolated states, the ground state 
and perhaps some isomers. In order to describe motion, 
one needs to consider a time-dependent mean-field the- 
ory, as e.g. time-dependent Hartree-Fock (TDHF). Large 
amplitude collective motion is related to low energy ex- 
citations, thus slow motion. This justifies the adiabatic 
limit known as adiabatic TDHF (ATDHF) [^3, E3| - It 
yields at the end, an optimized constraint Q and sub- 
sequent collective path {l^)} where q stands for con- 
tinuous series of deformations. The fully self-consistent 
optimization of the path is very cumbersome. It is a 
widely used approximation to use a simple quadrupole 
constraint as in (|Alb(l . The anomalies at large distance 
are avoided by a cut-off function / cut for which we use a 
Woods-Saxon shape 0|. The states are labeled with 
the dimensionless quadrupole moment (|AlcJ) which is 
rescaled with the total particle number A and the r.m.s. 
radius r. The index m can run over —2, —1, 0, 1, and 2. 
The path will be computed only along axially symmetric 
shapes corresponding to m = 0. 

The numerical solution is done by standard methods. 
Wavefunctions and fields are represented on an axially 
symmetric grid in coordinate space. An accelerated gra- 
dient method is used to iterate the single particle states 
(p n into their stable solution ^3] while the BCS+LN 
equations for v n are solved in each iteration step. An ex- 
tra iterative loop is included to maintain a wanted value 
of a 2 o 48] . 

Knowing the path, yields immediately the raw collec- 
tive potential as 



V(a 20 ) = E SHF (|$ Q20 )) 



(A2) 



11 



where -E S hf is the total SHF+BCS+LN energy for the 
given mean field state | ( E > Q20 ) includin g a cm. correc- 
tion as appropriate for the given force [23. The actual 
computations exploit explicit expressions in terms of the 
single-particle states (p„ and their occupation amplitudes 



2. Computation of masses 



The collective path allows to define a collective momen- 
tum operator as the generator of deformation P a |$ Q20 ) = 
id a20 \<S> a20 ) . At the same time, P a can also be interpreted 
as the momentum operator associated with collective dy- 
namics. The collective path {|$ Q20 )} is complemented by 
the dynamical response of the system by adding a dynam- 
ical constraint — /iP a to the mean-field equations yielding 
eventually a dynamical collective path. The adiabatic ap- 
proximation allows to handle the dynamical part in the 
linear regime, i.e. 



20Pq 



1 + iPa<3L C 



(dyn 



1$, 



(A3) 



The solution of the linear response equation thus ob- 
tained provides the dynamical response generator Q dyn ' > ■ 
Note that the whole energy functional is involved in 
the response. This is called self-consistent, or ATDHF, 
cranking, see e.g. |33l. l49|. The inverse collective mass 
for quadrupole motion is then obtained in straightfor- 
ward manner as 



B = - 



ld 2 E SHF (\<S> a20P J) 
2 



dp a " 



1 



= 2^ol[Qi dyn M^Qi dyn) ]]l^o)- (A4) 

The second form with the double commutator is not 
strictly applicable in connection with energy function- 
als (where the full H is not given) . It serves here only as 
a notational abbreviation to establish contact with stan- 
dard formulae for cranking masses. 

The same procedure is applied to the dynamical re- 
sponse to rotations. The collective momentum is here 
already known as it is the angular momentum, e.g J x . 
Solving the equations for the corresponding dynamical 
response yields the momentum of inertia as 



1 



id 2 E sliF (\$ a2 ^)) 

2 
1 

2 



dto 2 

= " ^ol^M^Q^]]!^), (A5) 



S( d y n ) 



carries 

(dyn) 



and similarly for y and z. The operator Qj 

the dynamical response in the same manner as Qa 
does that for the quadrupole motion. In practice, we are 
considering only axially symmetric shapes 0:20 ■ F° r then, 
we obtain 



For the case rotation, we simplify the response problem 
by computing the response with the stationary mean field 
h only (Inglis cranking). The approximation works very 
well for the considered cases. The critical region of small 
deformations does not contribute due to the topo- 
logical switching (|AT|) . 

For the full five-dimensional quantum corrections (see 
next subsection IA 3(1 . we also need to compute the inverse 
collective mass i? 7 and width A 7 for vibrations in the 7- 
direction. We do that for the vicinity of axial shapes 
by means of linear response. And we employ here the 
Inglis approximation using only the mean-field Hamilto- 
nian h in the response equations. All together, we have 
then the necessary ingredients concerning masses: the in- 
verse masses B, B 1 (— B), and the momentum of inertia 
B(a 20 )- 



3. Quantum corrections 

The GCM ansatz for the collectively correlated state is 
written as a coherent superposition over the path. How- 
ever, the states of the path correspond to wave packets 
in quadrupole space rather than to eigenstates of (X20 
(and similarly for the dynamical extensions in p a and 
^crank)- Thus they contain spurious contributions from 
collective motion which contribute to any expectation 
value. The strongest effects are found in the energy ex- 
pectation values which constitute the raw collective po- 
tential i|A2|l . These spurious contributions need to be 
subtracted. That is what one calls the quantum correc- 
tions or zero-point energies (ZPE) [T^. 

The correction for spurious center-of-mass motion is 
already part of the standard SHF scheme. The most 
important for the collective dynamics is the vibrational- 
rotational correction. These need to be considered as one 
entity because vibrations and rotations are closely con- 
nected pieces of the nuclear quadrupole topology. The 
recipe for purely axial vibration and rotation was given 
in [l2l f3T| . A recent model calculation has confirmed 
that ansatz and proven that the correction provides also 
a very good approximation to angular-momentum pro- 
jection, again for heavy nuclei. Here we want to account 
for the whole five-dimensional quadrupole dynamics (see 
next section FA 51 . Thus we are using the properly gener- 
alized rotational- vibrational correction 
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where J x y is the angular momentum about the x- or y- 
axis. The widths are the same for x and y because we 
evaluate everything at axial symmetry, i.e. at the point 
7 = 0. The total ZPE is decomposed as 



Vo) = (5-4. 9 (^)k:l ( A7a ) 

'A rotAF (ZPE) 



9(a) 



5ldxa{x 2 -l)e a ^-V 
J^dxe^ 2 - 1 ) 



(A7b) 



There is also a correction from spurious particle number 
fluctuations. This is already taken into account in an ap- 
proximate manner by the Lipkin-Nogami scheme added 
on top of the BCS pairing. 

All together, the quantum-corrected collective poten- 
tial reads 



V(a 2 o)=V(a 20 )-E^ E) (a 20 ) 



(A8) 



This is the quantity entering the collective Hamiltonian. 
The masses are associated with the collective kinetic en- 
ergies which are already of second order in the collective 
momenta. The quantum corrections on masses would 
correspond to terms of fourth order and are neglected. 



4. Retuning the particle number 

All states along the collective path are tuned to have 
the same average proton and neutron number. The ener- 
gies are corrected by approximate particle-number pro- 
jection at the level of the LN scheme. But the BCS 
states from which the collective path is composed still 
carry these particle-number fluctuations. As a conse- 
quence, the coherent superposition of the states along the 
path may change the average particle number again. One 
needs to readjust the correct average at the level of the 
collective dynamics |5l). To that end, one builds the col- 
lective picture of the particle number operator N — N in 
precisely the same manner as it was done for the Hamilto- 
nian. One obtains a particle-number potential, particle- 
number masses for quadrupole as well as triaxial motion, 



and particle-number contributions to the inertia. The 
expressions are the same as above with H replaced by 
N. The collective image of N is added as a constraint in 
the collective Schrodinger equation. 



5. The collective Schrodinger equation 

Axially symmetric quadrupole deformations are la- 
beled by a2o- The full space of quadrupole deforma- 
tions is explored when considering all a 2m with m <E 
{—2,-1,0,1,2}. This is convenient for spherical vibra- 
tor nuclei as it implies automatically the correct num- 
ber of vibrational degrees of freedom. It is, however, 
not well suited for deformed nuclei because rotations 
look rather involved in that frame. It is customary to 
transform by appropriate rotation into an intrinsic frame 
where <X2±i = and °22 = a 2-2- This defines three Eu- 
ler angles i9 as rotational coordinates. The remaining 
two relevant deformation coordinates 020 and 022 are ex- 
pressed in terms of total deformation /3 and triaxiality 
7 as 020 = /3cos(7) and 022 = /?sin(7)/V2 [53, 
Each triaxiality 7 which is an integer multiple of 60° 
corresponds to an axially symmetric shape. The cases 
7 = 0°, 120°, and 240° correspond to prolate axial defor- 
mations while 7 = 60°, 180°, and 300° are oblate. Rel- 
evant information is contained in one 60° sector of the 
plane, e.g. in the segment 7 G [0°,60°]. The other seg- 
ments can be reconstructed by axis exchange of principle 
axes. This symmetry under axis transformation has im- 
portant consequences for the representation of wavefunc- 
tions and potentials in the collective Schrodinger equa- 
tion: One has to obey mirror symmetry under 7 — > —7 
and axis-rotation symmetry under 7 — > 7 + 120°. The 
five-dimensional volume element d 5 a reads in the /3-j- 
frame d 5 a = (3 A \ sm(?>^)\ df} d'y d 3 9 where the 9 are the 
three Euler angles for the transformation from the lab- 
oratory frame into the intrinsic frame. In the following 
we will use the notation d 5 a as shorthand for the lengthy 
right hand side. 

The collective Hamiltonian has the form of a Bohr- 
Hamiltonian DjJS] generalized to /3-7-dependent masses 



fl-(coll) 

1 _ 3 
e fc (/3, 7 ) " 4 sin 2 ( 7 -f k ) 



-l-dfsB(l3 n )l3%- 1 
p 4 p l sin 37 



9 7 i? 7 (/3, 7) sin 37 <9 7 + 



L' k 2 



k=l 



20 fc (/3, 7 ) 



(A9a) 



X((3)+X(-/3) t X(f3)-X(-P) 



cos(3 7 ) Xe{B,V} 
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fee {1,2,3} (A9c) 
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(q\LM) = * iM (/3, 7 ,^) = J2^W)xi M (7,0<) f = 0,l,. 
{X° (7,^),^ = 0,l,...} - 0{c0B(3iAy)D^} = 



^^(cos(3 7) ) 



W M (7,^)^ = 0,1,...} = cos(A 7 )^ -(-l) A%3 sin(A 7 ) 



D (2) , D (2) 



(AlOa) 
(AlOb) 
(AlOc) 



A = 3 



+ v%2 + 1 
I 



where O means ortho-normalization of the set, [...] the 
integer part of a fraction, and a%b the modulo of a with 
respect to b. The operator L' k denotes the angular mo- 
mentum in the intrinsic frame. The deformed SHF cal- 
culations provide input along axially symmetric defor- 
mations a2o < 0- The collective dynamics needs to be 
performed properly in all five quadrupole degrees of free- 
dom. No strong peaks or wells in the 7 direction are to 
be expected for the nearly spherical or weakly deformed 
soft vibrator nuclei which we will consider here. It is 
thus an acceptable approach to interpolate the axial mi- 
croscopic results into the full (3-j plane. While potential, 
inverse masses and particle-number masses could be in- 
terpolated in a straightforward manner by l|A9bf) . the 
moment of inertia has three components. Respecting the 
collective symmetries we reconstruct them from the axial 
information as l|A9c|) . The deduction from axial infor- 
mation implies that we neglect the (3-"f coupling in the 
kinetic energy, i.e. Bp- — 0. The collective particle- 
number operator 7?( co11 ) is composed in the same form 
with the corresponding particle-number masses, poten- 
tials and moments of inertia inserted. 

The collective Schrodinger equations reads now 



(All) 



where the correction of the Fermi-energy ep o11 is to be 
adjusted such that J d 5 a ^+A>( co11 )* = N . 

The dynamics is formulated in the whole /3- 7 -plane but 
all necessary information is contained already in a seg- 
ment of 60°, as discussed in section lA~5l In order to meet 
the inherent symmetry conditions [lj, the wavefunctions 
of a + and a 2 + state are expanded in a symmetrized 
base IjAlOf) where the base mode v — determines the 
overall (3 dependence and the higher v shape the profile in 



the 7 direction. The Hamiltonian is very soft in 7 such 
that few j/-terms suffice for convergence (two or three, 
never more than 5). The D { $ K (f) are the well known 
Wigner D-Functions describing the rotation matrices for 
a state with angular momentum L |53|. It is notewor- 
thy that the structure of the rotational energy is that for 
the most general case where the considered nuclei have 
no special symmetry. For that reason the 2 + state must 
be a sum over all possible z-projections of the angular 
momentum which are K = 0, ±2 in the intrinsic frame. 

The remaining collective equation for the components 
i/j^((3) is solved numerically with standard methods. The 
wavefunctions and fields are represented on an equidis- 
tant grid in /3. Gradient iteration is used to find the few 
lowest eigenvalues and states. 



6. Computation of observables 

The solution of the collective Schrodinger equation, as 
outlined in the previous section, provides the energies 
directly. Expectation values and transition moments of 
other observables O need yet to be computed. The steps 
are, in principle, the same as done before for the energy 
(= Hamiltonian H respectively) and for the particle num- 
ber N. One has first to determine the collective image 
of the given observable O — ► O co ii and computes then 
the expectation value of that image with the collective 
wavefunctions [12. Actually, we do that for the compu- 
tation of the transition probabilities + — ► 2 + , the pro- 
ton B(E2) values. The observable here is O = Q20,prot- 
In case of the isotope shifts the observable is the collec- 
tive mapping of the radius. At present, we approximate 
the collective image by the raw expectation value only 
and neglect the kinetic corrections. 
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